First-principles prediction of high Curie temperature for ferromagnetic bcc-Co and 
bcc-FeCo alloys and its relevance to tunneling magnetoresistance 
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We determine from first-principles the Curie temperature Tc for bulk Co in the hep, fee, bec, 
and tetragonalized bet phases, for FeCo alloys, and for bec and bet Fe. For bcc-Co, Tc=1420 K is 
predicted. This would be the highest Curie temperature among the Co phases, suggesting that bec- 
Co/MgO/bcc-Co tunnel junctions offer high magnetoresistance ratios even at room temperature. 
The Curie temperatures are calculated by mapping ab initio results to a Heisenberg model, which 
is solved by a Monte Carlo method. 
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In the past few years we are witnessing a compelling 
race of different research groups^*^^ hunting the maxi- 
mum tunneling magneto-resistance ratio (TMR) of mag- 
netic tunnel junctions^ (MTJs) at room temperature 
(RT). MTJs, made of two ferromagnetic electrodes sepa- 
rated by an insulating barrier, open vistas to a wide field 
of technological applications, in particular in non- volatile 
magnetic random access memory, or new type of record- 
ing heads for ultrahigh-density hard-disk drives. Large 
temperature stability and other attributes are required 
which can be met by a high TMR at RT. 

After the first reported observations of TMR at room 
temperature^ in A^Oa-based junctions, reaching 30% 
at 4.2 K and 18% at 300 K, the achievement of high 
TMR values developed into a challenging materials 
science problem. Guided by theoretical predictions of 
extremely high TMR over 1000% for Fe/MgO/Fel, giant 
values of over 150% at RT have been experimentally 
achieved for fully epitaxial Fe(100)/MgO(100)/Fe(100)2 
and Fc(100)/MgO/Co(100)^ as well as poly-crystalline 
FeCo/MgO/FeCoi and FeCoB/MgO/FeCoB(001) 
MTJs4 Theory predicts even higher TMR values (at 
temperature T = 0) for epitaxial Co/MgO/Co junctions- 
than for Fe/MgO/Fe. Bcc-Co is a metastable phase and 
cannot be grown as a single crystal. Recently, however, 
it has been possible to grow bec Co in contact with 
MgO£±&±i and epitaxial bcc-Co/MgO/bcc-Co(100) 
MTJs show a high TMR, unusually stable^ with T. 

In most MTJs there is a large difference of the TMR 
between cryogenic and room temperature. The TMR, 
i.e., the relative change of resistance when the two fer- 
romagnetic leads are coupled either ferromagnetically or 
antifcrromagnetically, depends on the details of the elec- 
tronic structure, such as the spin polarization of the fcr- 
romagnet at the Fermi level Epr^ or the spin-dependent 
symmetry of the states at Ep^^- At T > 0, magnetic 
excitations mix the two spin channels, decreasing the 
TMR^ One can infer that a high Tc and spin stiffness 
entail a temperature stability of the magnetic structure 
and of the TMR. In this paper we investigate the Curie 
temperature of bcc-Co and Fei-^Co^ alloys and relate 
it to fcc-Co, hep-Co, and bcc-Fe, concluding that a high 
(calculated)^ Tc of bcc-Co is responsible for the temper- 



ature stability^ of TMR in bcc-Co/MgO/bcc-Co junc- 
tions. 

We calculate the Curie temperature using a standard 
recipe: the adiabatic approximation for the calculation 
of magnon spectrai 15 ! 16 Ab initio total-energy results are 
mapped to the classical Heisenberg model, 



H 



(1) 



Here, Jy are the exchange constants between the mag- 
netic moments at sites i and j, and &; is a unit vector 
along the moment of atom i. Tc was calculated within 
this model by a Monte Carlo method (using 1728 atoms 
in the supercell) by locating the susceptibility peak. 

The ab initio results are calculated within the gen- 
eralized gradient approximation^ to density-functional 
theory. We employ the full-potential Korringa-Kohn- 
Rostoker (KKR) Green function method^, using the co- 
herent potential approximation for the electronic struc- 
ture of the Fei-^Co^ random alloys. The exchange con- 
stants Jij are related to the Green function via the Licht- 
enstcin formula^ within perturbation theory assuming 
an infinitesimal direction change of the magnetic mo- 
ments at sites i and j. For Fe and Co we also compare 
to a complementary method, where finite-angle static 
magnon dispersion relations, J(q), are calculated on a 
dense mesh of (f-points in the Brillouin zone within the 
full-potential linearized augmented plane-wave method 
(FLAPW)i^ From J(q), the real-space constants Jy are 
found via a Fourier transform.— The two methods give 
the same trends. We stress that the trends are im- 
portant for our conclusions, and not the absolute val- 
ues of Tc, which are off the experimental value (bec Fe: 
T c = 1043 K; fee Co: T c = 1403 K)^l by -10%. 

Our calculations include bee Fe and Co, tetragonal- 
ized (bet) Fe and Co (considering growth on MgO(OOl)), 
the disordered alloys Feo.75Coo.25, Feo.50Coo.50, and 
Feo.75Coo.25, and the ordered alloy FeCo (in the CsCl 
structure). In order to see the effect of compression on 
the magnetic properties, we calculated the results at two 
lattice parameters: a=2.30 au=2.804 A, corresponding to 
the calculated equilibrium lattice parameter of bec Co, 
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FIG. 1: (color online) Density of states (solid line) and mean- 
field ordering temperature (dashed line) for bcc Co (left) and 
bcc Fe (right). 



and a=2.40 au=2.857 A. We find for hep Co an equilib- 
rium lattice parameter a = 2.488 A, with c/a = 1.623 
(experimental values are a = 2.51 A, c/a = 1.623), 
for fee Co a = 3.52 A (exp. value is a = 3.54 A), 
for bcc Co a = 2.804 A (extrapolated exp. value is 
a = 2.819 A)^ 2 . and for bcc Fe a=2.825 A (exp. value is 
a=2.867A). For Co and Fe on MgO(OOl) we find a mis- 
match of 5% in the surface lattice parameter (compared 
to -\/2aMgO = 5.63 A), resulting in a ratio c/a=0.857 for 
bet Co and c/a=0.909 for bet Fe. The ground state spin 
moments per Co atom, M hcp = 1.6, M fcc = 1.65, and 
M bcc — 1.75 /is, change with the crystal structure by 
less than 10%. The density of states (DOS) of bcc Co, 
shown in Fig. Q] exhibits the typical bimodal behavior of 
the bcc lattice with bonding and antibonding states. The 
spin polarization at Ep, P = —82%, is negative and has 
a larger (absolute) value than for bcc-Fe (P = 52%). Our 
calculated band structure agrees with the one of [23[ . 

Our results on the magnetic properties are summarized 
in Table |TJ A striking effect is that the exchange con- 
stants and the Tc increase with the Co concentration in 
Fei-^Co^- alloys, with a maximum of about Tc=1670 K 
for the ordered FeCo alloy, and a Tc ~ 1400 K for bcc 
and bet Co. In contrast, bet Fe suffers a 20% decrease of 
Tc compared to bcc Fe. We also find interesting trends in 
the magnetic moments and in the lattice-parameter de- 
pendence of the exchange constants. As the Co concen- 
tration increases, the local moment of Fe becomes larger, 
climbing from Mp ~ 2.2 (1b in pure Fe to about 2.6 [Ib 
in Feo.25Coo.75. The Co moment is comparatively inde- 
pendent of concentration, about 1.73-1.80 \ib- Moreover, 
the Fe-Fe first-neighbor exchange constants J\ (Fe-Fe) are 
strongly dependent on the lattice parameter for low Co 
concentrations x, but much less so for high x; the Co- 
Co and Fe-Co exchange constants, J\ (Co-Co) and Ji (Fe- 
Co), are much less affected by the lattice parameter at 
any x. The same behavior is found for the longer-distance 
exchange constants (not shown here). This leads to a 
strong dependence of Tc on the lattice parameter for Fe- 
rich alloys, and a weak dependence for Co-rich alloys. 

The change of the Fe properties upon increasing the 
Co concentration can be explained by observing that Fe 
progressively changes character from a weak fcrromagnct 
(i.e., both majority- and minority-spin d-states, S and 



cr , are only partly occupied for pure Fe) to a strong fer- 
romagnet (eft states become fully occupied); Co itself is 
a strong ferromagnet. The progressive change occurs be- 
cause the d(Fe)-d(Co) hybridization is weaker than the 
d(Fe)-d(Fe) hybridization, as the Co d states are more 
localized (they are deeper in the atomic potential well) 
than the Fe d states. As a result, an increased coordi- 
nation with Co results in an increase of the Fe moment 
and exchange splitting, driving the Fe & states lower and 
the Fe d* states higher in energy. The Co moment, on 
the other hand, is rather stable (around 1.8 /xb), since 
the Co d) states are always fully occupied^ The average 
moment per unit cell, M, first rises upon alloying with 
Co, due to the rapid increase of the Fe local moment, 
then peaks and drops for higher Co content (because of 
the comparatively lower Co moment). This peak of M in 
Fe-Co alloys is long know n 25 ' 26 ' 27 ' 28 , and is related to the 
triangular shape of the Slater-Pauling curve at the center 
of the 3d series 2 ^; density-functional calculations repro- 
duce this behavior, although they underestimate the local 
Fe moment by about 10% in these alloys^ 

In a weak ferromagnet the magnetic properties can 
be easily affected by structural perturbations (compres- 
sion or tetragonalization) because an extended part of 
the Fermi surface has d) and d^ character, allowing for 
d^-d^ charge transfer upon the structural perturbation. 
Such d^-d^- transfer can affect the local moment, but 
even more the spin susceptibility (and exchange con- 
stants Jij) which is sensitive to d'-d^ Fermi-surface cross- 
ing or nesting 2 ^ because of virtual tf-d^ spin- flip excita- 
tions. In the case of strong ferromagnetism (Co or FeCo 
alloys), small structural perturbations cannot cause a d^- 
d^ charge transfer, because the d) band is fully occupied 
and well under Ep. Thus the exchange interactions are 
more stable in Co and FeCo alloys. 

We proceed to the discussion of the exchange inter- 
actions and Tc by introducing^ the coefficient Jo = 
y^yn Joi 1 corresponding to the band-energy cost for flip- 
ping the magnetic moment of a single atom, reflecting 
a "single-site spin stiffness". It is related to the mean- 
field Curie temperature via ksT^ F = 2Jo/3 (fee is Boltz- 
mann's constant). (It is well-known that mean-field the- 
ory overestimates Tc, but it is a useful tool for trends 
analysis.) By treating Ep as a parameter, we calculate 
T^ F (E) as a function of band filling; in this way we are 
able to see the individual contribution of the states at 
each energy to the exchange interactions. In Fig. [1] we 
show T(f F (E), together with the DOS, for bcc Fe and 
Co. Evidently, T^ lF (E) has a very similar form for Fe 
and Co, up to the approximately rigid band shift (it also 
has slightly higher values for Fe because of the stronger 
d-d hybridization); thus, Fe and Co can be compared in 
a unified picture. The negative values of T^ F (E), peak- 
ing around —2 eV for Co and —1 eV for Fe, indicate 
an antiferromagnetic coupling, known from (5-Mnj2! At 
higher energies the double-exchange mechanism sets in 
and T(f F (E) obtains strong positive contributions as E 
crosses the final part of the d^ states and the anti-bonding 
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Alloy 


a(A) (c/a) 


M Fo 


Mco 


M 


Af oxp [25j J^- tc 


7U0-U0 


jrFc — (Jo 






Fe bcc 


2.804 


2.107 


— 


2.107 


2.22 


17.45 


— 


— 


970 


1120 




2.857 


2.197 


— 


2.197 




15.25 


— 


— 


900 




Fe bet 


2.978 (0.909) 2.321 


— 


2.321 





12.71 


— 


— 


750 




Feo.75Coo.25 


2.804 


2.421 


1.785 


2.263 


2.45 


25.80 


24.54 


31.78 


1390 




(disordered bcc) 


2.857 


2.512 


1.803 


2.335 




26.93 


24.61 


31.67 


1490 




Feo.50Coo.50 


2.804 


2.523 


1.778 


2.151 


2.35 


28.55 


22.05 


30.76 


1600 




(disordered bcc) 


2.857 


2.588 


1.801 


2.194 




27.99 


22.13 


30.27 


1600 




FeCo 


2.804 


2.729 


1.727 


2.228 


2.42 




- 


28.84 


1660 




(ordered CsCl struct.) 


2.857 


2.803 


1.743 


2.273 






- 


28.50 


1670 




Feo.25Coo.75 


2.804 


2.560 


1.759 


1.959 




26.64 


20.01 


28.33 


1520 




(disordered bcc) 


2.857 


2.633 


1.791 


2.002 




26.42 


20.06 


28.01 


1540 




Co bcc 


2.804 




1.751 


1.751 






18.48 




1420 


1670 




2.857 




1.790 


1.790 






18.71 




1370 




Co bet 


2.978 (0.857) 




1.740 


1.740 






16.76 




1380 




Co fee 


3.519 




1.646 


1.646 






13.82 




1280 


1200 


Co hep 


2.487 (1.623) 




1.602 


1.602 






14.66 




1300 


1350 



TABLE I: Calculated (within KKR) magnetic moments (in Hb), first-neighbor exchange constants Ji (in meV), and Tc(in K) 
of bcc FeCo alloys, bcc Fe and Co, bet Fe and Co, fee and hep Co. The bcc structures are considered at two lattice constants, to 
show the effect of compression. For the bet structures, the in-plane lattice constant was adapted to MgO and c/a was relaxed. 
Tc is calculated by a Monte Carlo method with exchange constants derived from the KKR or FLAPW methods. 



d^- states. Finally it drops to zero once the d^ states are 
filled and the exchange mechanism is no more present. 

In Fe, Ep is located at a steep, ascending point of 
Tq(E), before the maximum. Therefore, small struc- 
tural perturbations resulting in band shifts have a strong 
influence the T c of Fe, as seen in TableQ] For Co, T^ F (E) 
is already descending at Ep, but is not as steep as for 
Fe. From this argument, Co is expected to have a higher 
and more robust Tq than Fe. In a FeCo alloy Ep is in- 
between, at the maximum of T^ F (E); then one expects 
the highest and most robust Tc, as is found by the Monte 
Carlo calculations. 

In conclusion, ferromagnetism is found to be more ro- 
bust in bcc Co than in bcc Fe. The calculated Xc(Co) 
is higher and stable with respect to structural changes, 



even in the tetragonalized (bet) structures. This can 
be advantageous in the temperature dependent TMR 
of Co/MgO/Co MTJs compared to Fe/MgO/Fe, as ob- 
served in recent experiments*^ The Tq of bcc Co is cal- 
culated to be the highest among all Co phases. Although 
the highest Tc was found for FeCo alloys, they have the 
drawback that the TMR can drop if FeCo disordered 
layers appear at the interface with MgO, because the 
symmetry, enforcing spin polarization of the current, will 
be absent. Taking the excellent coherent spin-dependent 
transport properties in conjunction with the MgO(100) 
barrier into consideration, Co/MgO/Co may be an ideal 
junction for high TMR values at elevated temperatures. 

We gratefully acknowledge enlightening discussions 
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